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Abstract 



The temperature dependence of the anomalous sector of the effective action 
of fermions coupled to external gauge and pseudo-scalar fields is computed 
at leading order in an expansion in the number of Lorentz indices in two 
and four dimensions. The calculation preserves chiral symmetry and confirms 
that a temperature dependence is compatible with axial anomaly saturation. 
The result checks soft-pions theorems at zero temperature as well as recent 
results in the literature for the pionic decay amplitude into static photons in 
the chirally symmetric phase. The case of chiral fermions is also considered. 
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I. INTRODUCTION 



It is well established in the literature that the axial anomaly has a temperature indepen- 
dent form [p]-[T3||. This result is consistent with our present understanding of anomalies, 
since they are induced by the ultraviolet divergences present in the theory whereas the finite 
temperature modifies the infrared sector only, namely, by imposing periodic or antiperi- 
odic boundary conditions in the Euclidean time direction 0-16]. In particular for Weyl 



fermions, the known topological origin of the anomaly guarantees its independence under 
changes on the space-time manifold Ul7 |. Likewise, the axial anomaly is mass independent 
and also density independent [[TJ|. Analogous statements hold for the parity anomaly in 
odd-dimensional theories fl8|| . 

Recently, it has been found in |B| that the anomalous amplitudes, such as the neutral 
pion decay into two photons, are temperature dependent. At first, this would seem surprising 
since the numerical value of the pion width is very well accounted for by the axial anomaly 
prediction and furthermore, such decay is considered the standard proof that chiral anomalies 
are not just mathematical artifacts. In order to clarify this point it should be noted that 
even at zero temperature the axial anomaly and the neutral pion decay are different entities; 



Adler's theorem pOfl , relates the pion decay to the axial anomaly in the soft-pion limit only. 
Translated to the language of effective actions, this means that, at leading order in a gradient 
expansion, the pseudo-parity odd component of the effective action [^1|] (or equivalently, 



the component that is imaginary in Euclidean time or also the component containing the 
Levi-Civita pseudo-tensor) is accounted for by the gauged Wess-Zumino-Witten (WZW) 
action |22]j2T[| . In other words, all other contributions are of higher order. They contain 
more gradients and are therefore sub dominant in the soft-pion limit. Such higher order 
terms will necessarily be chirally invariant since the gauged WZW action saturates the 
anomaly. If one tries to write down Euclidean (Lorentz) and chiral invariant terms of the 
same order as the WZW action, it is immediately clear that they vanish identically. At finite 
temperature the conditions are less restrictive and the pion decay is no longer determined 
by the axial anomaly [^,@]. Indeed, at finite temperature the time direction is privileged 
and Euclidean invariance is partially broken. Moreover, the effective action is not expected 



to admit a gradient expansion with local terms in general [ 16,23-25 1 . This lower symmetry 



allows to write down new chiral invariant terms which can compete with those coming from 
the anomaly rendering the anomalous (or better, pseudo-parity odd) amplitudes temperature 
dependent even at leading order. 

Ref. ]19| makes use of a linear sigma model with constituent quarks to study the mesonic 
decay amplitudes at temperatures near the chirally symmetric phase through direct com- 
putation of the relevant Feynman diagrams. As already noted, a non trivial temperature 
dependence is found and furthermore the neutral pion decay turns out to be suppressed in 
the chirally symmetric phase. This conclusion is confirmed in |26| using functional methods. 
Both calculations use the imaginary time formalism to introduce the finite temperature. A 
similar calculation is carried out in |27| using the real time formalism. There it is found that 
the pion decay amplitude is indeed temperature dependent although the suppression in the 
chirally symmetric phase is not reproduced. In [^] the problem is studied in full generality 
regarding the kinematical conditions of the pion and photons within the real time formula- 
tion. The analysis there indicates that the discrepancy comes from the different kinematical 
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configurations assumed, namely, static photons in [19| versus on-shell photons in ||27|| . It 
is also noteworthy that, according to [28|], the vanishing of the pionic decay amplitude in 
the chirally symmetric phase will presumably be recovered when quantum pionic fluctua- 
tions are properly taken into account since they regulate the infrared sector. Another set 
of calculations of anomalous (and non anomalous) mesonic amplitudes should be mentioned 
here [^[23"H 25|j . They correspond to the low temperature limit and thus well within the 
phase were chiral symmetry is broken. In this case the spirit of chiral perturbation theory 
applies. These calculations make use of effective Lagrangians of pions. In this approach, the 
temperature dependence comes through pionic loop corrections (instead of quark loops as 
in the constituent quark models cited above) and the pionic Lagrangian itself is assumed to 
be temperature independent. 

In the present work we carry out a calculation of the pseudo-parity odd component of 
the effective action of fermions in two and four dimensions in the presence of external non 
Abelian vector, axial and pseudo-scalar fields on the chiral circle at finite temperature, and 
at leading order in a suitable expansion. We will use the ^-function prescription and will 
emphasize the chiral symmetry preservation aspects such as the anomaly and explicit vector 
gauge invariance. Some implications for chiral fermions are also presented. 



II. GENERAL CONSIDERATIONS 

The Euclidean action describing fermions in the presence of external bosonic fields is 
/ 'ipDifj, where D is the Dirac operator 

D =1) + 4 75 + MW" , (1) 

— 0^ + V/j, is the covariant derivative, U 15 stand for U in the subspace 75 = +1 and U^ 1 
when 75 = —1. M is a constant c- number mass term, the constituent mass of the quarks, 
Vn(x), A^ix) and U(x) are matrices in flavor (and color) space. M is real and positive, V^(x) 
and Afj,(x) are anti-Hermitian, and U(x) is unitary. The matrices 7 M and 75 are Hermitian 
and satisfy 

{7m> 7^} = 2<V, , 75 = Vdlo ■ ■ ■ ld-i , (2) 

d = 2,4 being the space-time dimension and 772 = h — 1- We will use the imaginary 
time formalism to implement the finite temperature condition, namely, through periodic 
boundary conditions in the Euclidean time direction for the bosonic fields and antiperiodic 



conditions for the fermionic fields with period (3 = 1/T, T being the temperature [14-16 
We will assume a space-time topology of the form R x S 1 . As a consequence Euclidean 
(Lorentz) invariance is partially broken due to the anisotropic boundary conditions. 
Chiral symmetry corresponds to the transformation 

d - n^n^Dng*n? > (3) 

where Pr,l = |(1 ± 75) are the projectors on the subspaces 75 = ±1 and ^h,l(^) are 
independent unitary flavor matrices. In terms of the chiral fields V^ L (x) = ± A M , the 
Dirac operator can be written as 
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D = {If R +MU)P R + [fl L +MU- l )P L , (4) 

where D R,L = + V R,L are the chiral covariant derivatives, and the chiral transformation 
takes the form 

Df L - n R \D*> L n R , L , u^ni l un R . (5) 

Vector gauge transformations correspond to the diagonal subgroup Q R — Ql, i.e. D — > 
f2 _1 Dfi. For finite rotations, Q r< l — ^ x p{<^r,l) the matrices a R ^ being ant i- Hermit ian. 
Infinitesimal vector and axial rotations are defined by 5a R ,L = Say ± 5a a- In the Abelian 
case finite axial transformations form a group. 

The fermionic effective action is formally defined as W = — Tr log(D). To renormalize 
the ultraviolet divergent trace we will adopt the ^-function prescription pO|-0^1, that is, 
W = — J2 n 1°S ^n\s=o- Here, A n are the eigenvalues of D and s = is understood as an 
analytical continuation from sufficiently negative s. This method is mathematically well 
founded |3(J and it is particularly convenient in the context of finite temperature since it 
automatically preserves vector gauge invariance under both small and large gauge trans- 
formations [|33],[n|. This is because, within the ^-function prescription, the effective action 
depends only on the spectrum of the Dirac operator which is invariant under gauge trans- 
formations. Note, however, that in general the effective action is invariant modulo 2ni, due 
to the multivaluation of the logarithm. 

It is convenient to introduce the concept of pseudo-parity transformation [ 21f1 , namely, 



— > Vp, M — > M, — > — and U — > U~ l . Due to parity invariance, which involves 
and additional (x , x) — > (x , — x) in the fields, the pseudo-parity odd component of the 
effective action is that containing the Levi-Civita pseudo-tensor. Also, it corresponds to the 
imaginary part of the effective action. As it is well-known, the real part of the effective action 
may be regularized preserving chiral symmetry and the chiral anomaly is only essential (i.e., 
not removable by counterterms) in the pseudo-parity odd component |T7| . Here we will 
concentrate on this latter component, which will be denoted W~ . 

The consistent chiral anomaly, defined as the variation of the effective action under an 
infinitesimal chiral rotation, is given by 



5Wj =2 = % - J tr (F - A 2 ) 5a A 



SW. 



d = 4 ' 12vr 



l — J tr (3F 2 + F\- AAFA - {F, A 2 } - A 4 ) Sa A , (6) 



where we have adopted a standard differential geometry notation: D = D^dx^, V = V^dx^, 
A = A^dXfn, F = D 2 = dV + V 2 , Fa = {D, A}, the dx^ anticommute and dx^dx\ ■ ■ ■ dxd-i = 
£m...d-id d x. It is noteworthy that the anomaly depends on the gauge fields only and not on M 
or U(x) and also that there is no anomaly associated to purely vector gauge transformations. 

The anomaly can be integrated to yield the gauged WZW action P3|,ZI,|53-3B[|, which in 
two dimensions takes the form 

r(V, A, U) = / tr(i? 3 ) + / tr (v R R - V L L + V R U~ X V L \J - V R V L ) , (7) 

where R = U~ 1 dU and L = UdU^ 1 . In the Wess-Zumino term, / tr(i? 3 ), the integral 
refers to a three-dimensional manifold on the gauge group which interpolates between the 
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original U(x) configuration and a fixed configuration belonging to the same homotopy class. 
The latter can be taken as a constant at zero temperature. The finite temperature case 
is discussed in great detail in Ref. |J7|. There it is shown, in particular, that T does not 
depend on small variations (in the topological sense) of the interpolating path and also that 
the variation of F is independent of the homotopy class representative chosen. For this 
invariance it is essential that the Wess-Zumino integrand is a closed form. 

Alternatively, the same gauged WZW action can also be written as [|T^,[^-fID],|23] 

r(V, A,V) = ~Jtt Qi? c 3 - R C F R - 2AF R + ^ - p.p.c , (8) 

where R c = R — V R + U~ 1 V L U, F R = D\ = dV R + V R transform covariantly under Q R 
and are invariant under fl L , and p.p.c stands for pseudo-parity conjugate, i.e., A — ► —A, 
U — > U^ 1 (and thus exchanging the labels R with L everywhere). This version can be 
obtained directly by starting from the Wess-Zumino action and applying minimal coupling, 
R — > R c . The result is then chiral invariant but it is no longer a closed 3-form, i.e., it has 
a spurious dependence on the particular interpolating path taken. This is cured by adding 
new terms involving the gauge fields and preserving vector gauge invariance. An advantage 
of this version is that the chiral breaking terms are manifestly polynomial thus yielding a 
polynomial anomaly. These terms cannot be removed by local polynomial counterterms, 
however, since they do not form an exact three-form by themselves. The corresponding 
expression in four dimensions is 

r(V, A, U) = - J tr(^i? c 5 + 2R C F 2 R - R*F R + R C F R U^F L U + 2AF R F L 

+AAF R - 8A 3 F R + ^-A 5 ) - p.p.c . (9) 



(The LR versions are given in section f7T[ ) 

We remark that the mathematical statement is that the infinitesimal variation of r(V, A, U) 
equals the consistent anomaly given in the right-hand side of eqs. (|6]), provided that V, A 
and U transform under chiral rotations as in eq. @. Therefore, U can actually stand for 
any field taking values on the gauge group and transforming as f^f/fi/j. This is true when 
U is the pseudo-scalar field appearing in the Dirac operator, but would apply also for a U 
suitably constructed out of V and A, for instance. This remark is specially relevant when 
M — 0, since then there is no pseudo-scalar field in the Dirac operator. In any case, by 
construction, the gauged WZW term r(V, A, U) and the effective action W~(V, A, U; M, T) 
differ by chiral covariant terms only. 

In order to carry out the ulterior calculation we will fix the chiral gauge: by taking 
Q R = Q and Ql = UQ a new field configuration is obtained such that U = 1. Let us denote 
the new gauge fields by V M and A^. The field Q(x) represents a vector gauge freedom. It will 
be convenient to partially fix this gauge by imposing stationarity of the Vo field. As shown 
in Ref. [I8|1, the remaining vector gauge freedom consists of two kinds of transformations, 



namely, stationary gauge transformations and discrete transformations of the form £7 = 
exp(a;oA(a;)), where A(x) is any anti-Hermitian matrix commuting with Vq(x) and with 
eigenvalues of the form 2nin/ (3, for integer n. The quantization of A ensures the preservation 
of the periodic boundary conditions. In the Abelian case, the discrete transformations are 
large in the topological sense. 
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The difference between the effective actions of the original and rotated configurations 
is accounted for by the gauged WZW term, which vanishes identically when U(x) = 1, 
therefore 

W~(V,A,U;M,T) = W~(V, A, 1; M, T) + T(V, A, U) . (10) 

By construction, since the chiral gauge has been fixed, the effective action of the (V, A,U = 
1) configuration defines a chiral invariant action when re-expressed in terms of the original 
fields. In the actual calculation, which necessarily truncates the exact result to some order, 
chiral invariance will be maintained provided that the calculation in terms of the rotated 
fields preserves vector gauge invariance. 



III. CALCULATION OF THE EFFECTIVE ACTION 

In this section we will explain our calculation of the effective action at finite temperature 
and the expansion used. The result of the calculation is discussed in the next section. Let 
us introduce the ^-function 

n s (D)=Tr(D s )=Y,K, (11) 

n 

which is ultraviolet finite if Re(s) < —d and a meromorphic function of s with simple poles 
at s — —d, —d + 1, . . . , — 1 |50|] . Then, the ^-function prescription is 

W(D) = ~n.(D) . (12) 



ds 



s=0 



Applying Cauchy theorem 



n s = -Ti I ( 13 ) 

h 2ni D - z K J 



where the integration path in the complex plane starts at infinity following for instance the 
negative real axis, encircles the origin clockwise and goes back to infinity along the same 
ray p0|j40|] . In this regularization, the mass-like term z is responsible for the explicit breaking 



of chiral symmetry needed to allow for the chiral anomaly. 

To deal with the trace on the space-time degrees of freedom, we apply the Wigner 



transformation method [4C.1S1: 



n s = - [ f^-W I — ^tr<0l— — ^ 10>. (14) 

Here, p = ik, p = u n with cu n = 2iri(n + |)//3 and the sum on n refers to all integers, 
tr refers to flavor and Dirac spaces. |0) is the state of zero momentum and zero energy in 
the space-time Hilbert space, normalized as (sc|0) = 1. In particular, this implies c^|0) = 
(0|c^ = whenever a derivative inside D reaches any of the ends of the matrix element. 
Also, (0|0) = Jd d x. 

In order to proceed we will work with a chiral gauge fixed configuration, i.e., with the 
Dirac operator 
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D = ^(D ll + A tll5 ) + M, d V = 0. (15) 

(For simplicity, in this section we will not use a special notation for the rotated fields. 
Also, for the correct interpretation of the formulas, we remark that in this section we are 
not using a differential geometry notation.) The main idea of the calculation, in order 
to be able to carry out the sums and integrations indicated in eq. fll4|) , is to perform a 
series expansion on some of the pieces contained in the Dirac operator. Depending on the 
precise choice made one can obtain perturbation theory, gradient expansions, inverse mass 
expansions, etc. We want to organize the calculation so that operators of higher dimension 
are also of higher order, however, it is essential to do this preserving gauge invariance. In 
the Wigner transformation formulation, the integral over k projects out the contribution 
which is invariant under time-independent gauge transformations whereas the sum over 
frequencies does the same thing for the discrete transformations |40|JT8|] . (Loosely speaking, 



a time-independent gauge transformation can be compensated by a shift in k and a discrete 
gauge transformation can be compensated by a shift in u n .) Therefore, gauge invariance 
is preserved by expanding in powers of ~fD+ 4 75 while keeping j> +70-D0 + M — z in the 
denominator. At finite temperature, making a further expansion in powers of D would break 
gauge invariance because D transforms discretely under discrete gauge transformations. 
The deep reason for this is that counting powers of do effectively means to study the change 
of the functional under dilatations in the time coordinate, and such dilatations are not 
supported by the boundary conditions. 

Treating D fully non perturbatively in the way just described is possible since, due 
to our choice of chiral gauge fixing, the quantities p M , d , V and M — z appearing in the 
denominator are all commuting. (See [0 for such a non-perturbative treatment in three 
dimensions.) However, in the spirit of retaining only lower dimensional operators, a further 
natural expansion is that in powers of the quantity Dq defined by Dq(X) = [Do, X] which 
preserves gauges invariance. 

Let us illustrate this kind of calculations by detailing the procedure in the two dimen- 
sional case. Of course, due to the identity 7^75 = ie^v, in this case it is algebraically 
simpler to use a single non anti-Hermitian vector field — — ie^A^ in the calculation. 
Nevertheless, for greater similarity with the four dimensional case, we will keep the vector 
and axial fields as independent variables . Since the leading order terms are in principle of 
the form VqA and A V, we have to retain terms of first and second order in the expansion 
of the ^-function, 



n.,i*» = -f~T,f ^tr<o| f - 1 (jD+ 4 75) 

J 2tt (3 ^ Jr 2m \ >y Q + jp + ji 



loQ + IP + A* 



1 ( lD +4 l5 )——± — —( lD +4 l5 )——± — — ]|o>, (16) 



loQ + JP + H loQ + jp + fi 70Q + IP + V 

where Q = u n + D and /1 = M — z. Since Q, p and /1 commute with each other we can use 
the identity 

1 fi - 70Q - IP 



loQ + 1P + V 



A = /2 2 -Q 2 + k 2 . (17) 



The Dirac trace can then be evaluated. Retaining only the pseudo-parity odd terms one 
finds 
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r dk 1 r d7 1 / /I 1 

0- l+2 = 2z y -- E / r — ^ s /^tr(0| — ( + QA + (3fc 2 - /i 2 ) (a -D + £>-A 

-Q^,i[D, Q] + [D, Q}^A Q + A Q^DQ + QD^QA^j I|0) . (18) 

Next, all the "naked" .D, i.e., those which are not inside commutators, are brought to 
the right of the matrix element. For instance 

^l = i D + i {g ' [Ag]} l- (19) 

This produces commutator terms of the form [15, X] which are gauge invariant under time- 
independent gauge transformations. There are also non-covariant terms, of the form tr(0|X£)|0), 
which can be replaced by tr(0|XV|0), since <9j|0) = 0. Thus at the end all appearing oper- 
ators are multiplicative in a?-space, that is, x becomes just a parameter and the operators 
effectively act on a Hilbert space which is the tensor product of flavor and time spaces only. 

The problem now is that the sum over the variables k, u n and z cannot be done in a 
straightforward manner because they appear in different operators which do not commute. 
At this point, one can insist on obtaining a strict expansion in powers of D and A^ (i.e., 
without further expanding in powers of Dq) ||18[| . Instead of doing so, we will simplify the 
problem by retaining only terms which are of lowest dimension in fields and derivatives. The 
counting is defined as follows: M will be taken of order zero, V^, A^ and count as first 
order each. This is equivalent to count the number of Lorentz indices. The treatment of D 
requires some care: the quantity Q is in principle of zeroth order due to u n , however, inside 
commutators u n does not contribute, so Q becomes D and counts as first order. Then 
(recalling that expanding in powers of D is forbidden) the strategy to follow is to bring 
all naked Q (the explicit ones as well as those inside A) together to the left of the matrix 
element. There they commute and the indicated sums and integrals can be easily carried 
out. In doing this we keep only leading order terms. Moving explicit the Q operators to the 
left can be done systematically, generating commutator terms. On the other hand, moving 
A operators to the left is more subtle in the sense that it cannot be done exactly in closed 
form. It can only be done up to some given order in the expansion. This comes from the 
identity 

X± = ±X-2%[D ,X\± + I [Do, [D , (20) 

which shows that new terms with a A" 1 factor at the right are generated. Nevertheless, such 
terms always come with commutators and thus they are of higher order than the original 
one. Therefore repeating the commutation operation a sufficient number of times, all naked 
Q will end up at the left, modulo higher order terms. Once this has been achieved, all Q 
not at the left has been replaced by D and they appear inside commutators. The next 
thing to observe is that in all the naked Q (explicit and implicit) which are now at the left, 
D can be replaced by Vo since (0|<9o vanishes and Vq is stationary. So finally we have a 
multiplicative operator both in x- and xo-spaces and hence (0|X|0) is just / d d xX. Note 
that there is no difficulty of principle in computing the expansion at any given order. In 
practice, it is simpler to bring the naked Q to the left and the naked D to the right doing 
the truncation at the same time. 
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Applying this method to eq. (|18D , the pseudo-parity odd ^-function in two dimensions 
becomes, at leading order, 



^ r dk 1 ^ r dz . 

([D,A ]+2A„V 




A 2 



10), (21) 



where effectively Q = u n + Vq. In this expression, the contributions from [D, A ] and A V 
vanish after the integration over k. The cancellation of non-covariant terms of the form 
tr(0|XV|0) is a non trivial check of the calculation. Likewise, the contribution from [A, Do] 
vanishes using integration by parts and cyclic property. (The cyclic property does not apply 
directly here since the operation (X) = £) n tr(0|X|0) is not a trace, however, it can be shown 
that (XY) = (YX) whenever X is multiplicative in a>space and Y is a function of D [[TH]-) 

The only remaining term QA~ 2 A presents a further subtlety. In principle, this term 
would be of order one but it is actually of second order. This is because, after summing 
over frequencies, the factor QA~ 2 yields an odd function of Vo (recall that do is no longer 
present here) and thus this factor is of first order. In general, the Levi-Civita pseudo-tensor 
requires the saturation of d Lorentz indices in a (^-dimensional space-time, thus the leading 
order term will be of order d. 

In four dimensions the calculation of the leading contribution is algebraically more in- 
volved, but it can be done along the same lines. For instance, there appear terms of the form 
(A~ 2 [[A, Dq\, F]), which in principle would count as fourth order, however, using again the 
cyclic property, it is clear that the non trivial contribution from A starts at second order 
in Vo, so this terms is actually of sixth order. (This was to be expected in this and similar 
formally ultraviolet divergent terms, since otherwise there would be an spurious contribution 
to the scale anomaly which is absent in the pseudo-parity odd sector.) 

The calculation can alternatively be done using a (d + 2)-dimensional formalism, which 
relates the pseudo-parity odd part of the effective action in d dimensions with the baryon 
number in d + 2 dimensions |il- fj5| , [i0| (for instance, the two dimensional Wess-Zumino 



term (127r) -1 / tri? 3 is the correctly normalized baryon number in four dimensions). Such 
formalism gives directly the effective action as the integral of a (d+ l)-form, as in eqs. (§,0). 
One advantage of this procedure is that enforcement of vector gauge invariance is sufficient 
to fix the action without introducing a ^-function regularization, but working in more space- 
time dimensions is also more involved. The results found with this formalism (at least in 
two dimensions) are consistent with those found in the (^-function prescription. 

Let us remark that fixing the chiral gauge by U — 1 and OoVq = is essential to carry 
out the calculation. Fixing Vo to be stationary allows to set do = in the operators Q 
moved to the left. Eliminating U allows to use the identity in eq. (|TT|). An exception is 
the case of Abelian and stationary fields, since in this case U, Vo and do commute and this 
allows to carry out the calculation without fixing U = 1. Such calculation can be done most 
conveniently using the (d + 2)-dimensional formalism and in particular it checks, once more, 
that the anomaly is temperature independent. 
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IV. THE EFFECTIVE ACTION AT FINITE TEMPERATURE 



Recalling the exact formula eq. flTU|), we define the leading order of the pseudo-parity 
odd component of the effective action as 

W^(y, A, U; M, T) = W£ d .(V, A, 1; M, T) + r(V, A, U) , (22) 

where V and A refer to the rotated fields. The chiral invariant term, W^^V, A, 1;M, T), 
has been computed with the method described in the previous section and is given by 



WW,c. = ~ 2i j tr Ki- 4 ) > ( 23 ) 
W^ cX = J tr faftiA, T) + + 8 - V ^) MA, J= A ) - f ^ 



for two and four dimensions, respectively. Here, a differential geometry notation has been 
used with A = Aidxi, J- = \[D h Vj\dxidxj and J- a = \Ph Aj]dxidxj, and Aq actually 
stands for Aodx . Further, the coefficients are functions of Vq (and M and T) given by 



<n id) 



a '( 2r /(^5?/ r ^ (M -^)„„' (24) 



where Q = uj n + Vo (with uj n = 27rz(n+ §)//?) and A = (M — z) 2 — Q 2 + fc 2 . Besides, a factor 
<ixo is implicit when m is odd so that in this case (pfb is a 1-form. When m is even, * s 
an even function of Vo and is of zeroth order, on the other hand for odd m, (fffa is odd in 
Vo and it counts as first order. As a consequence all the terms in W^ lead are of dimension d. 
The following simpler formulas are equivalent whenever they are convergent: 

^ - fJ—LiT^ (25) 





are 




r dk 




1 2^ 




r dk 


"i 


' 2tt 



dXr 



These functions are not all independent. The following relations are useful 

¥>?J = - + MVS ■ (27) 
So the four-dimensional action can also be written as 
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WWi = / tr ( 2 ^{ A ^ + (J M Vg! - 7^) A{ A J~ a\ + y (~^| - MVS) -4 

(28) 

In two dimensions, c } contains only a term with the structure VqA.. The absence 

of AqV can be understood as a two dimensional peculiarity. In this case the Dirac operator 
depends only on the two combinations Vo — iA\ and Vi + iAo and not on the four fields 
independently. Because gauge invariance requires D\ + iAo to appear only inside commu- 
tators, the possible terms of dimension two of the form AqV would have to come from 
[T>i + iAo,Vi + iAo] which is identically zero. On the other hand, in four dimensions it 
is not clear why W^, &d ci contains no terms with the structure AqVVV, and why no time 
derivative should appear. 

The gauge invariance of W{^ &d c i is obvious under time-independent gauge transforma- 
tions. Under discrete gauge transformations Vo transforms to Vo + A, so its eigenvalues are 
shifted by an integer multiple of 2m/ '(5. Due to the sum over frequencies, are periodic 
functions of Vo with period 2iri/ (3 and the result is invariant. 

Another important test is the zero temperature limit of the result. As noted before, soft- 
pion theorems imply that at zero temperature, the only pieces present at leading order in 
the pseudo-parity odd sector are those coming from the gauged WZW term, which already 
saturates the anomaly. Indeed, the new contributions would be chirally invariant, and 
thus gauge invariant in terms of the rotated fields (i.e., in the chiral gauge U = 1). The 
possible leading order terms consistent with Euclidean and gauge invariance are {J 7 a) m two 
dimensions and (J 7 J 7 a) and (J 7 aA 2 ) in four dimensions, all of which vanish identically. Of 
course, there are new chiral invariant terms at zero temperature at sub-leading orders, for 
instance, in two dimensions, the fourth order term is (e.g. using the formulas in Ref. [j40| ) 



J (fxe^Mi^ - [A^A v })[D x ,Ax]) . (29) 



127T M 2 

Therefore, restoration of Euclidean invariance requires W^ &A ; to vanish in the zero tem- 
perature limit. In this limit / ^y.-i J2 n becomes / ?fJ^- All terms containing with 
odd m vanish: since u n is now a continuous variable, it can be shifted to eliminate Vo- The 
resulting integral vanish due to antisymmetry of the integrand. This cancellation is illus- 
trated by eq. ( p6|) . Likewise, the coefficient of the term Aq{A,J-a\ also vanishes at zero 
temperature. This can be seen using the form in eq. (|28|) and computing ip^l with eq. (pop. 



V. ANOMALOUS AMPLITUDES 

In this section we will obtain the mesonic amplitudes derived from the effective action 
at finite temperature. We will consider only the Abelian and stationary case, since it allows 



to compare with previous calculations in the literature Hl9] , Pq| . 

In the Abelian and stationary case, V = V and A = A — hd(p, where U = exp(0), thus 
the chiral invariant part of the action at leading order reads 

u -,,,„.,„ = -2» / y?S (a - , (30) 

WWc, = / (*P$F + (y MVS - ^) AoF A ) (a - \d<P 



11 



In order to find the amplitudes, let us retain the leading order in an expansion in powers of 
Vo- A simple calculation, gives 

^l = -^(l-W) + C>(V 3 ), (31) 



where we have introduced the dimensionless function 

/=E « ss , (32, 

which has limits / — > 1 when T — > and / — >• when M — > 0. On the other hand, in the 
Abelian and stationary case the gauged WZW term reduces to 

T 2 (V,A,U) = ^J V dtf> 

r 4 (V, A,U) = --^[ (3V F + A F A ) dct> . (33) 
Therefore the leading anomalous amplitudes can be read from the following effective actions 

W2 = hl v " + 2(1 ~ f)A) ' (34) 

^ J (3V F + A F A ) (/d0 + 2(1 - . 



W 
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These formulas suggest that the temperature dependence of the amplitudes are independent 
of the space-time dimension. In the two cases considered, introducing a finite temperature 
amounts to make the substitution d<p — > fdcf) + 2(1 — f)A in the anomalous amplitudes at 
zero temperature. In particular, the term VoFd<f), related to the process tt° — ► 77, gets an 
extra factor of / as compared with the zero temperature amplitude. For M<T, this factor 
behaves as 

7 M 2 ( M A \ 

> = 4?« ! V +o (t0- (35) 

A result which is in agreement with those in Refs. Hl9jp6|. 

In the zero temperature limit / = 1 and the amplitudes reduce to the gauged WZW 
term, as they should. Thus in this case there is a tight relation between the anomalous 
amplitudes and the anomaly. On the other hand, at finite temperature the amplitudes are 
modified, yet the anomaly is preserved since the terms added to the gauged WZW action 
are chirally invariant. 

The fact that / vanishes as M — > implies that the amplitudes involving <p cancel in 
the chirally symmetric phase. In principle, such a cancellation is to be expected on general 
grounds (and thus it is a non trivial test of the calculation): since the variable U no longer 
appears in the Dirac operator when M = 0, the effective action should also be independent 
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of U. More generally, the effective action should depend analytically on the external fields 
V/j,, A^ and the combination MU. There are, however, two considerations which affect this 
conclusion. First, the argument may be spoiled by infrared divergences. For instance, at zero 
temperature the dominant anomalous term is the WZW action which depends on U rather 
than MU. At finite temperature, u n is always different from zero so no infrared divergences 
should appear in the fermionic effective action. Second, even if the argument holds for 
the exact effective action, it needs not apply when only the leading order is retained. The 
truncation might introduce a spurious <fi dependence at M — 0. The situation at M = 
beyond the Abelian and stationary restrictions is further studied in the next section. 

Another important remark is that the fact that / — > as M — > does not directly 
imply that the amplitude for Tc-f-f vanish in the chirally symmetric phase. For instance, 
in a linear sigma model, the relevant piece of the action is gip(M + a + i^^TTf)^ hence 
the pion field tt(x) is related to <fi as it ~ M0, as a consequence the relevant amplitude 
goes as f /MVoFdir and the renormalization factor due to the temperature is f/M. In our 
calculation, in agreement with [^9| , p6| , / = 0(M 2 /T 2 ) and f/M — > as M — > 0. However, 



in Refs. [p7| , p8| , / = 0(M/T). Whereas / still vanishes as M — > 0, the 7177 amplitude does 
not vanish, presumably due to the different kinematical conditions assumed [p8] . 



VI. MASSLESS FERMIONS 

In the case of massless fermions, the gauge fixing condition U — 1 is meaningless because 
there is no field U in the Dirac operator and no chiral rotation is necessary. Thus the 
calculation of the effective action described in section |T| can be applied to any gauge field 
configuration ( V, A) provided only that the requirement d V = is satisfied. In order not 
to unnecessarily complicate the notation, we will use the symbols V and A to denote, not 
the original fields, but the fields after a vector transformation so that Vo is stationary. That 
is, the Dirac operator is 

D = 7m (£> m + A m 7 5 ), d V = 0, (36) 

and the effective action at leading order is given by W 1 ~ ad _ c i in eqs. (p3|) but using V and A 
instead of V and A and with M = 0. The required massless functions can be computed 
in closed form and are given by 

^ = ~V 0i ^ 4 J = -ttVo, (37) 

Z7T on 

where V actually stands for the 1-form V dx . (A detailed analysis shows that in this 
formula V is to be understood modulo 27rz//5 so that periodicity is preserved.) Note that 
these functions no longer vanish as T — > in the massless case; the two limits M — > and 
T — ► do not commute. The leading order actions become 

= i /tr(VoA) , (38) 

7T J 

WW y > A ) = -7^2 / tr ( 3 MA F} + A {A, F A } + 2V A 3 ) . 
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Consistently with dimensional arguments, these actions are temperature independent (ex- 
cept, of course, through the boundary conditions). It should be noted that they are not 
really local polynomial actions due to the gauge fixing: when expressed in terms of the 
original fields, these actions are in fact non-local. 

In the massive case, the correct chiral transformation of the action was not an issue, since 
the chiral gauge had been fixed. In the present case, however, the correct transformation is 
not evident by construction. If these formulas are applied, not directly to (V, A) but to a 
chirally rotated configuration (V, A) (with Vo stationary) chiral symmetry would require 

W{^(Y, A) = W^XV, A) + r(V, A, U) , (39) 

where U denotes the rotation. This equality holds for the exact effective action, what is not 
evident is that it should hold for its leading order too. Note that an equivalent statement 
is that the right-hand side is independent of U, that is, there is no spurious dependence 
on U. The discussion of the previous section shows that this is true in the Abelian and 
stationary case. (In fact, the eqs. ( [31] ) which give ip^ and <^4 1 at lowest order in V , turn 
out to coincide with the exact ones when M — 0, eqs. (|37|) .) For Abelian but time dependent 
configurations, still A = A — \d(j), where U = exp(0), and V = V (there is no vector gauge 
transformation involved since 8qVq = <9oVo = 0). Then, in two dimensions, the right-hand 
side of eq. ( |3"9"D yields 

W^(Y,A) + ±f Vdo<j> (40) 

(where do stands for dodxo). This formula implies that in general there is a spurious depen- 
dence on U . 

In order to consider the general non-Abelian case, note that eq. ( |59"D is equivalent to 
say that W^ d {V, A) displays the correct chiral anomaly. A detailed calculation shows that 
these actions do not reproduce the correct chiral anomaly in general, but they do so for 
stationary configurations. It is rather remarkable that, even if restricted to the stationary 
case, these actions are able to saturate the non-Abelian anomaly in two and four dimensions. 
This is more so since the calculation of W{^ &d , detailed in section [III], knows nothing of the 
anomaly or the WZW term. From this point of view, it constitutes a non trivial check of 
the calculation. 

The fact that the anomaly is not always reproduced means that, in the massless case, 
the truncation of the effective action to its leading order violates chiral symmetry. This 
is because we have done an expansion in powers of the operator 7-D+ A 75 which is not 
chiral covariant; under chiral rotations it mixes with the other piece of the Dirac operator, 
7o-Do- This implies that the correct anomaly is only recovered through cancellations among 
the axial variation of terms of different order. (A different matter is directly computing the 
axial anomaly use our expansion. The axial anomaly within the ^-function renormalization 
is given by Tr(— 2a J 47 5 D s ) s=0 [fi4"l,fl0| which, in the massless case, can computed with the 
same technique described in section |TTT| and gives the correct result in closed form.) 

It is interesting to compare our result in the two-dimensional case with the exact effective 
action for the massless Dirac operator (the Weyl determinant) known in closed form 



The exact result follows from the observation that, upon complex analytical extension of 
the chiral group parameters, chiral transformations are sufficient to bring any gauge field 
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configuration to a space-time constant configuration. Therefore, the exact effective action 
is given by two terms. First, the gauged WZW term associated to the required analyti- 
cally extended chiral rotation and second, the effective action corresponding to the constant 
configuration. This latter term will depend on the temperature, as will also depend higher 
orders in our expansion. In the Abelian case (and concentrating from now on on the temper- 
ature independent term only) the corresponding exact effective action for the pseudo-parity 
odd sector is 

W-(V,A) = - ( V t A t . (41) 

71 J 

The labels £ and t stand for longitudinal and transverse, respectively, that is, 

V t = d{i^)- 1 d v V v ) , V t = V-V t . (42) 

This action is non-local thus it cannot be removed by a local polynomial counterterm, i.e. 
by a suitable choice of the renormalization prescription. In the stationary case, Vt = Vq 
and Ai = A, so the exact effective action coincides with that given in eq. (0). That is, in 
this particular case, the leading order is exact (in the temperature independent sector of the 
effective action). On the other hand, for non stationary configurations, an expansion of the 
exact result to extract its leading order is not well-defined due to the presence of d~ 2 in the 
definition of Vg, i- e -> due to infrared divergences in the expansion. 

The exact Weyl determinant is not known in four dimensions. The trick of analytical 
extension of the chiral group only doubles the dimension of the chiral group and thus it is 
insufficient to rotate to zero all the components of the gauge configuration. Nevertheless, 
in the Abelian case it is easy to write down a pseudo-parity odd action which saturates the 
correct anomaly, namely, 

W~(V, A) = J (2FAV t - FVA £ - F A AA e ) . (43) 

(A version of it exists for any even number of dimensions.) Unlike the two dimensional 
case, this action does not reduce directly to that given in eq. fl3~8|) for stationary configura- 
tions. This does not necessarily imply that they are inconsistent with each other, since the 
difference are terms which are subject to infrared ambiguities if one insists on a gradient 
expansion. 

Of course, for massless fermions it is much more natural to work with chiral fields. Since 
right and left fields are completely decoupled, the effective action must satisfy 

W~(V, A) = W~(V R ) - W~{V L ) + P(V R , V L ) , (44) 

where P is a local polynomial introduced by the renormalization prescription. In terms of 
the chiral field v (say V R ) the consistent anomalies take the well-known form 

SWz = ^- j tr (F - v 2 ) 5a , 
= — / tr vdSa , 

47T J 

5W4 = [ tr (-2F 2 + Fv 2 + vFv + v 2 F - v A ) 5a , 

Aon 2 J V ' 
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' r tr (Fv - -v 3 ) dSa . (45) 



24tt 2 J V 2 

Where, D = d + v,F = D 2 = dv + v 2 and 5t> = [Z?, 5a]. The corresponding gauged WZW 
actions are 

r^^^-^/tr^+^/tr^) 

= / tr (X 3 + t; 3 - 3(i? c + v)F) , 

r 4 («, E7) = -^5^ / tr (^ 5 ) + ^ / tr - WvR + v 2 R 2 - -vRvR + vR 3 ^j 

= ~ 24(U / tr + ^ " 5(jR ' + w3)F + 10(jRc + V)F2 ) • (46) 

Here, R = U~ 1 dU and R c = R — v, and Z7 transforms as UQ when Z) transforms as Q~ 1 DQ. 
In the (d + l)-dimensional version of formulas, the polynomial term is the normalized d+1 
Chern-Simons term and thus it yields the correct anomaly fl7| . 



The Abelian actions in Eqs. (|4]) and ( fHj|) satisfy the decoupling condition, eq. (0), with 
W - {v) = -i-J VtV , P=±JVA, 

W-(v) = ^Jv,vdv, P = ±- 2 j A VdV, (47) 

in two and four dimensions respectively. In each case the polynomial P is uniquely deter- 
mined by vector gauge invariance. These actions are constructed so that they saturate the 
anomaly, that is, they correspond to the gauged WZW term associated to the chiral rotation 
which brings the field v to the gauge vg = 0. (This is achieved by taking U = exp(</>) where 
4> is defined by vg = d(j>.) Of course, besides the gauged WZW term the exact effective action 
has another contribution (namely, the exact effective action of the purely transverse field) 
but it is chiral invariant by construction and, at least in two dimensions, it can be shown to 
be pseudo-parity even. The same procedure can be followed in the V A version of the theory, 
that is, taking the gauged WZW term associated to axially rotate to the gauge Ag = 0. In 
two dimensions this procedure reproduces eq. ( f4"T|) and thus it gives nothing new, however, 
in four dimensions it yields 

W-(V,A) = --±- ( (3FV + F A A)Ag. (48) 

1Z7T Z J 

This action differs from that in eq. (^) by a chiral invariant term of the form (FV t A t ). 
Since such term is not a local polynomial (and as far as I can see, is not identically zero) 
both actions are not related by a change in the renormalization prescription. In principle 
only the action in eq. ( |i"3|) , which satisfies the decoupling formula (44[), could come from 
integration of the fermions. 

The chiral version of our leading order action corresponds to take V R = v and V L = 
in the V A version, eqs. (j38|). This gives 

WW(«) = ^/tr(i*ti), (49) 
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which refer to the gauge 8qVo = 0. 

In the two dimensional case there is an important difference with the V A form of the 
action: due to the gauge condition, d v = 0, the only allowed infinitesimal gauge transfor- 
mations are time-independent ones. Under this restriction it can be immediately verified 
that the chiral two dimensional action at leading order always yields the correct anomaly; 
the configuration no longer has to be stationary as in the V A form. (In the V A form of 
the theory, the analogous statement holds if the formula is applied only in the chiral gauge 
d V = d A = 0. Eq. (|40D confirms this observation since in such a gauge <p can only be 
stationary and the spurious term Vd <f) cancels.) This means that for any configuration, 
the effective action can computed by integrating the anomaly up to a t> -stationary gauge 
and applying our leading order formula. The result will be independent of the particular 
^o-stationary gauge chosen. Unfortunately, the action so constructed differs from the cor- 
rect one by gauge invariant terms, unless both v and V\ are stationary. What would be 
needed is to "rotate" the original configuration to the "gauge" d v = d Vi = 0. This can 
be done using an analytical extension of the chiral group, in the same spirit as Rothe's 
calculation pE5f| . 

Not unexpectedly, in four dimensions the situation is worse. The anomaly is reproduced 
by W^lead. ( v ) if the configuration is stationary, but not otherwise. Therefore it is not possible 
to give a well-defined effective action in this case; the precise value would depend on the 
particular t> -stationary gauge taken. It would be necessary to bring the configuration to 
the stationary case, but analytical extension of the chiral group is not sufficient to do this. 

VII. SUMMARY AND CONCLUSIONS 

We have studied the temperature dependence of the effective action of fermions in the 
presence of external bosonic fields. Due to the complexity of the problem, only the leading 
order terms, that is, those which can compete with the anomalous ones, have been retained. 
However, the calculation could in principle be carried out to higher orders as well. The tem- 
perature dependence of the anomalous amplitudes is shown to be fully consistent with chiral 
symmetry including the known temperature independence of the axial anomaly. The known 
zero temperature limit, given by soft-pions theorems, is verified. Also, for temperatures near 
chiral symmetry restoration, previous results in the literature obtained with the imaginary 
time formalism for static final states are reproduced. Finally, some non trivial tests are also 
verified in the case of massless fermion, but full [/-independence of the effective action as 
M — > is only reproduced in the stationary case, due to the truncation at leading order. 
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